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Abstract 



We investigate the colhsional stabihty of magneticahy trapped uhracold 
molecules, taking into account the influence of magnetic fields. We compute 
elastic and spin-state-changing inelastic rate constants for collisions of the 
prototype molecule ^^©2 with a '^He buffer g function of the magnetic 

field and the translational collision energy. We find that spin-state-changing 
collisions are suppressed by Wigner's threshold laws as long as the asymp- 
totic Zeeman splitting between incident and final states does not exceed the 
height of the centrifugal barrier in the exit channel. In addition, we propose 
a useful one-parameter fitting formula that describes the threshold behavior 
of the inelastic rates as a function of the field and collision energy. Results 
show a semi-quantitative agreement of this formula with the full quantum 
calculations, and suggest useful applications also to different systems. As 
an example, we predict the low-energy rate constants relevant to evaporative 
cooling of molecular oxygen. 
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I. INTRODUCTION 



The probable success of experiments aimed at producing magnetically trapped ultracold 
molecular samples depends heavily on the effects of coUisional processes. For example, 
paramagnetic alkali dimers can be produced via photassociation of ultracold atoms , but 
the resulting molecules, typically in high-lying vibrational states, are subject to vibrational 
quenching collisions which can release a large amount of energy and dramatically affect 
the efficiency of the cooling. Alternatively, cold molecules in their vibrational ground states 
can be produced either by thermal contact with a cold helium buffer gas |^ or by Stark 
slowing, for species that possess an electric dipole moment Collisions are of obvious 

importance to buffer-gas cooling (BGC), as well as to forced evaporative cooling (EC) that 
will be required to lower the temperature of these gases further and achieve for instance 
Bose- Einstein condensation (BEC). Both processes require large elastic collision rates to 
thermalize the gas. 

So far EC has not been realized in practice for molecules, but the success of the BGC 
technique for the production of cold CaH |^ and PbO [0] molecules suggests that in the 
near future it will be possible to achieve BEC using cold molecules. This would open the 
way for a number of new and fascinating experiments. 

In order to be magnetically trapped, atoms or molecules must be in a weak-field seeking 
state, i.e. a state whose energy increases with the strength of the magnetic field. For 
each trappable weak-field seeking state there is in general a lower-energy untrapped strong- 
field seeking state, in which the molecules experience a force away from the center of the 
trap. Collisions can drive transitions between trapped and untrapped states. These "bad" 
collisions can cause heating or atom loss. It is therefore important to assess the rate constants 
for the inelastic collisions. 

In a series of previous papers the resilience of molecular oxygen against spin-changing 
collisions was investigated, in collisions of O2 molecules both with a helium buffer gas [§, 
and with other O2 molecules 0. These studies found that spin-changing rates due to 
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spin-rotation coupling could be quite large. However, in the case of the ^ molecule, 
where in the limit of zero field the only allowed exit channels are energetically degenerate 
with the incident channels, spin-flipping transitions require boosting the centrifugal angular 
momentum from L = to L = 2, meaning that these processes are strongly suppressed by 
the Wigner's threshold laws at collision energies smaller than the height of the exit channel 
centrifugal barrier. 

The results in Refs. [|,^ considered only the case of a vanishing external magnetic 
field, which is obviously not the case in experiments that trap molecules using spatially 
inhomogeneous magnetic fields. The present paper therefore explores the role that fields 
play in determining spin-changing collision rates. As we demonstrate, the presence of a 
magnetic field causes a Zeeman asymptotic splitting between incident and exit channels, 
thus lifting the collision energy higher relative to the centrifugal barrier in the exit channel, 
and removing the Wigner's suppression. 

Studies of spin-changing ultracold collisions in the presence of an external magnetic field 
have been performed so far only for atomic species |]TU|-[T^. In this paper, we present a 
detailed dynamical study at cold and ultracold temperatures for the atom-diatom system 
^^©2 — ^He in a field. The basic model is described in Sec. II. In Sec. Ill we calculate elastic 
and inelastic rate constants for collisions involving the lowest-lying trappable state of ^''02 
over a wide range of field values (from up to 5000 gauss), then discuss the dependence 
of the rates on collision energy for several representative values of the field. This system 
is of direct relevance to BGC of molecular oxygen. Generally, it allows us to quantify the 
removal of the Wigner's law suppression as the field increases in strength. On this basis we 
determine a simple one-parameter fitting formula that reproduces the trend with field and 
energy of the loss rates. In Sec. IV we use this formula to extend previous results on O2-O2 
collisions to estimate the influence of the field on EC of this system. 
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II. THEORY 



As mentioned in the Introduction, we will consider in this paper molecules consisting of 
two ^''O atoms, whose nuclear spin i is equal to 5/2. We assume that total spin I = ii + i2 is 
conserved in the collision and polarized to its maximum value 1 = 5, implying that the even 
molecular rotational states are separated from the odd ones |jTB|]. We limit the discussions 
of this paper to the "even-iV" manifold of molecular states, which is more appealing for 
cooling purposes [^,^ having a spin 1 paramagnetic ground state. 

The electronic spin S of the O2 molecule has magnitude S" = 1 in the electronic ground 
state we are concerned with throughout this paper. The angular momentum S is 
coupled to the molecular rotation angular momentum N to give J, the total molecular 
angular momentum, which assumes the values N — 1, N and + 1 for A^ > and is 1 for 
A^ = 0. The Hamiltonian operator H02 for molecular oxygen in the presence of an external 
magnetic field B can be written as: 

Ho, = 5eN2 + H/, + HB (1) 

where Bg is the rotational constant. The ^^02 molecule is considered to be a rigid rotor, 
with internuclear distance frozen to the equilibrium value of tq = 2.282 bohr (the rigid rotor 
model has been shown to be very accurate for this system at the investigated collision energies 
jl^). The fine-structure Hamiltonian H/s and the Hamiltonian for the interaction of 



the molecule with the external magnetic field follow the treatment in Ref. [|15|, disregarding 
the molecular hyperfine interaction. 



The fine structure Hamiltonian is given as |T5 

'2\i/2 



H^,= (^) AT2(S,S)-T2(a,a)+7N-S (2) 



where a is a unit vector parallel to the molecular axis and T is a second-rank tensor |16,17 



The fine structure parameters 7 and A have been taken from Ref. [0, where they have been 
determined by microwave spectroscopy. 
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The interaction of the field with the electronic spin can be expressed as [jT9|,E0|: 



Hb = 5- ■ /ioS ■ B (3) 

where B indicates the external magnetic field, g is the g-factor of the electron and /xq is 
the Bohr magneton. Following we ignore a small interaction between the field and the 
rotational angular momentum. 



The matrix elements for H/s and have been given in Ref. |jT5|, eqs. (A5) and (A6), for 
a Hund's case b basis set. We note here that the molecular rotational quantum number N is 
no longer strictly a good quantum number for the molecular states, because different values 
of are coupled together by the fine-structure Hamiltonian Hjg and by the interaction 
with the external field. The molecular total angular momentum quantum number J is still 
a good quantum number with respect to H/s, but not with respect to the field interaction 
term. However, its projection Mj on the quantization axis is still conserved. 

Consequently, our basis functions should be labeled as |nMj), where n is a shorthand 
index denoting the pair of quantum numbers {N, J) in the field dressed basis |T^. However, 
the coupling between different A^'s is weak (the fine-structure coupling is small compared 
to the rotational separation) and A^ can be considered "almost" a good quantum number. 
Similarly, J is also approximately good for laboratory-strength magnetic fields, so that we 
can use without confusion the label \NJMj). 

Magnetic trapping is strongly related to the behavior of the molecules in a magnetic 
field. The low-energy Zeeman levels of oxygen are displayed in Fig. 1 for the even-A^ 
species. (Throughout this paper we report energies in units of Kelvin by dividing by the 
Boltzmann constant fc^. These units are related to wave numbers via 1 K = 0.695 cm^^). In 
order to be trapped in the usual magnetic traps a molecule must be in a weak-field-seeking 
state, i.e., one whose energy rises with increasing magnetic field strength. Thus the state 
\NJMj) = |01 1) is the lowest-lying trappable state of the even-A^ manifold, and this is 
the state on which we focus our attention below. This state is indicated by a heavy line in 
Figure 1. Higher- lying states with N > 2 are energetically forbidden at low temperatures. 



It is clear from Fig. 1 that for any trapped state there is an untrapped, strong-field- 
seeking Zeeman state at a lower energy. These states are not merely untrapped but anti- 
trapped, experiencing a force away from the trapping region. In a magnetic trap collisions 
with buffer gas atoms, or more generally with other molecules, will therefore ultimately 
deplete the trap of its molecular population. The time available for cooling processes like 
BGC or EC, as well as the lifetime of a molecular Bose-Einstein condensate, is therefore 
limited and knowledge of the rate constants for spin-flipping collisions is essential to predict 
their feasibility. 

In Ref. PI the theoretical framework for atom - diatom scattering was derived in the 
limit of zero external field, along the lines of the model originally due to Arthurs and 



Dalgarno pT| , p2| , and properly modified to incorporate the electronic spin of the oxygen 
molecule. Here, the formulation of the scattering problem is further extended to account for 
the interaction with the external magnetic field. 

The full Hamiltonian operator describing the He - O2 collision is given by 
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2/i 



+ Uo, + V{R,e) (4) 



after multiplying the wave function by R in order to remove first derivatives. Here fi is the 
reduced mass for the He-02 system, R is the modulus of the Jacobi vector joining the atom 
to the molecule's center-of-mass, is the centrifugal angular momentum operator, and H02 
is the molecular oxygen Hamiltonian defined in The potential term V, depending on 
both the Jacobi vector R and on the bending angle 6 that the molecule's axis makes with 
respect to R, accounts for the He - O2 interaction. We use the ab — initio potential energy 



surface (PES) by Cybulski et al. p3[, which approximates the true well depth to ~ 20%. 

The full multichannel calculation requires casting V{R, 9) in an appropriate angular 
momentum basis. Our field dressed basis for close-coupling calculations is then 

\02{^T.-))\Yle{'S))\NJM.jLMLM) (5) 

where the electronic spin quantum number S is not explicitly indicated being always equal 
to 1 in the problem treated here. The quantum number L stands for the partial wave 
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representing the rotation of the molecule and the He atom about their common center 
of mass, Ml is the projection of L onto the laboratory axis, and M. is the laboratory 
projection of the total angular momentum, M. = Mj + Ml- At the collision energies of 
interest we assume that the oxygen electronic state and the helium atom state are preserved 
and therefore we suppress the first two kets in equation in the following. 

We note here that, at variance with the formulation in the zero field limit, the total 
angular momentum of the system J' (equal to N + S + L) is no longer a good quantum 
number, because different J7's are coupled together by the interaction with the external 
field. This means that the dynamical problem is no longer factorizable for different values 
of J', thus requiring larger numbers of channels to be treated simultaneously. The problem 
is still factorizable for A4, but in general the number of channels to be included for each 
calculation is much larger than in the previous case. Numerical details of the calculations 
will be given in the next section. 



The coupled-channel equations are then propagated using a log-derivative method ||2| 
and solved subject to scattering boundary conditions to yield a scattering matrix S: 

{NJMjLMl\S{M)\N'J'M'jL'M'^) (6) 

As already noted, the projection of the total angular momentum M. is still a good quantum 
number, implying that Mj + Ml = Mj + Mj^. Note that in general each of the quantum 
numbers A^, J, Mj, L, and Ml are subject to change in a collision, consistent with conserving 
M. 

Following Ref. P], the state-to-state cross sections are given by: 

aNjMj^N'j'M', = E \{NJMjLMl\S - I\N'J'M'jL'M'l)\^ (7) 

f^NJMj lMlL'M'^ 

and the corresponding rate coefficients are given by 

KnJMj^N'J'M'j = VnJMjC^NJMj^N'J'M'j (8) 

where vnjmj is the relative velocity of the collision partners before the collision. For no- 
tational convenience we will in the following refer to collisions that preserve the incident 



molecular quantum numbers as "elastic," and those that change the quantum numbers as 
"loss." 



III. RESULTS 

In this section we consider elastic and state-changing, inelastic rate constants for the 
incident channel \N J Mj) = |0 1 1). At the investigated collision energies (from 1 fjK up 
to 10 K) there are two open inelastic channels, namely J Mj) = |0 1 0) and |0 1 — 1), 
both of which are untrapped (see Fig. 1). 

A. Magnetic field dependence 

We begin by computing rate constants for the elastic and spin- flipping transitions |0 1 1) 
— > |0 1 — 1) and |0 1 1) — > |0 1 0) in the low-field hmit. Results in this section refer to 1 
HK collision energy and are converged using partial waves up to L = 6 and including the 
rotational states A'^ = 0, 2, 4 and 6 for the oxygen molecules. The maximum value of R to 
which the coupled-channel equations are propagated depends on the strength of the field, 
ranging from 600 bohr in the case of smallest field to 450 bohr for highest values. These 
parameters assure rate constants convergent within less then 5 %, which is adequate for our 
purposes. 

In principle, the scattering matrix should be determined for each possible value of the 
projection of the total angular momentum Ai. However, we know that s-wave collisions 
dominate the incident channel at ultralow collision energies, which corresponds for our in- 
cident channel to the value A4 — 1. We have verified that including in the calculation only 
the Ai = 1 channels changes the results by less than 1 % at /xK energies; in this section we 
therefore include only this contribution. The number of channels to be propagated according 
to the given convergent quantum numbers is then only 205. 

Fig. 2 shows elastic and inelastic rate constants at energy E = l/xK and for low values 
of the field. The inelastic rate constants are non zero even at zero field, as shown in Ref. 
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P]. However, in this limit the final states are degenerate with the initial state, and inelastic 
transitions are strongly suppressed by the presence of a d-wave centrifugal barrier in the exit 
channel, whose height is about 0.59 K. This effect is able to suppress the molecule loss, at 
least as long as the collision energy does not exceed the barrier height [Q. 

As soon as a field is applied, the thresholds are no longer degenerate in energy, so that 
the energy in the exit channel is not as far below the centrifugal barrier. As a consequence, 
inelastic transitions are no longer as strongly suppressed, even in the limit of very low 
collision energy. Rather, they increase dramatically even in a weak field, with rates being 
boosted by 5 or 6 orders of magnitude in a 1 gauss field. On the other hand, elastic scattering 
is nearly unchanged by the field. 

This sudden increase of the inelastic transition rates can be reproduced semiquanti- 
tatively by applying the distorted wave Born approximation (DWBA) p2[, as has been 
successfully done for the magnetic dipolar interaction of cold alkali atoms [p!0| , p!I| as well as 



in a number of problems in cold collisions p5|j26|] . The first order DWBA is a simple two 
state perturbation approximation applicable in cases where inelastic scattering is weak in 
comparison with elastic scattering. The DWBA expression for the off-diagonal K-matrix 
elements is: 

K,f = -n / MR)V,f{R)ff{R)dR (9) 

J 

where fi and ff represent the energy normalized scattering wave function of the initial and 
final states calculated on the diabatic potential corresponding to the states involved in the 
inelastic transition. V^/ is the corresponding diabatic coupling term of the Hamiltonian (^. 
The first-order DWBA result is shown as a dotted line in Fig. 2. 

The DWBA also yields information on the threshold dependence of the loss rates on 
energy and field. To see this, first note that the spin-changing processes we are considering 
are strongly dominated at low collision energy by s-waves in the incident channel, and by 
d-waves in the exit channel. This change of partial wave is necessary to conserve angular 
momentum during a collision that changes the molecule's spin. For small values of the mag- 
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netic field (for wliich tlie Zeeman splitting does not exceed the height of the exit channel 
centrifugal barrier) the exit channel is still in the threshold regime, whereby the wave func- 
tions fi and // in (^ can be approximated by the small-argument limit of energy-normalized 
spherical Bessel functions, 



/, cx ^JkijLXkiR) oc ihRf^^'/^ ff oc ^kfjL,{kfR) oc (fc/i?)^^+l/^ (10) 

where fcj and kf are the incident and final wave numbers and U and // are the incident and 
final partial waves. Assuming a small- i? cutoff to insure convergence of the integral in (^ 
with respect to singularity in the coupling potential Vif{R), the energy-dependence of 
the K-matrix element is 

oc kf^^'/'k'/^'/' (11) 

By considering the relationship between Kif and the effective rate constant K^jmj^n'j'm'j-i 
it is straightforward to show that the rate constant behaves approximately as 

Knjmj^n'j'm'^ « {E + AMjg^oBf'-^^'^ (12) 

where E is the collision energy, and we have taken into account that the final kinetic energy 
in the exit channel is incremented by an amount AEb = AMj g ^qB corresponding to the 
linear Zeeman shift. AMj (= Mj — M'j) stands for the difference between the initial and 
final values of Mj in the two channels involved in the transition. The actual value of AEb is 
modified by quadratic Zeeman shift, but the linear approximation is adequate for achieving 
a simple fitting formula. In the present case these shifts amount to less than 10% changes 
in the approximated rate constants. From eq. ([T2|), considering that in our case Li = and 
Lf = 2, a simple expression for the rate constants can be derived: 

fE + AMJg^^oB\'/' 

J^NJMj^N'J'M'j = J^O y J [i-o) 

where Kq represents an overall scaling constant and Eq is conveniently chosen as the height 
of the centrifugal barrier in the exit, d-wave channel. In the limit of very low collision energy, 
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the AMjQfioB term obviously dominates over E, leading to a nonzero rate constant, as is 
the case for exothermic collisions. 

This simple expression allows us to interpret the threshold behavior of the rates with the 
field, explaining the 5/2 exponential dependence on B found in our calculation and shown 
explicitly in the bilogarithmic plot for the rate constants (Fig. 3). Here the dashed lines 
represent a fit to the rate constants in the limit of zero magnetic field, yielding coefficients 
Kq = 2.73 X 10"^^ cm^ sec~^ and Kq = 1.45 x 10~^^ cm^ sec~^ for the transition to 
the final states |0 1 — 1) and |0 1 0), respectively. Apart from zeros in the actual rate 
constants, the overall trend is indeed Kj^jm,-^n'J'm' oc 5^/^. The zeros in the real rates 
arise from interferences between the s-and d-wave radial wave functions, as we have verified 
qualitatively by the DWBA. Nevertheless, the simple one-parameter expression (p!3D provides 
a reasonable upper bound to the complete calculation which, it will be recalled, requires a 
calculation involving 205 coupled channels. 

The simple formula (p!3| ) holds, of course, only when both incident and final channels are 
in the threshold regime. Assuming low incident energies, this restriction therefore limits the 
size of magnetic field for which (|TB|) applies. Namely, this expression is only useful when the 
Zeeman energy splitting between incident and final states remains smaller than the height 
of the d-wave centrifugal barrier. For the channels considered here, these fields are 2430 
gauss and 4860 gauss for the |01 — 1) and |010) final states, respectively; Vertical arrows in 
Fig. 3 indicate these field values. Relation (p!3D serves as a useful quick fitting formula for 
data at smaller fields, and allows us to generalize the results obtained in this paper also to 
different systems. Possible apphcations are discussed in section IV. 

When stronger values of the magnetic field are considered, the DWBA is no longer able 
to reproduce the full calculation, because the coupling is no longer a weak perturbation. 
Strong field interaction mixes up different channels leading to a much more complicated 
picture which can not be explained in terms of a simple two state model. However, in this 
limit the loss rate constants seem to be unacceptably large anyway, owing to the effectiveness 
of spin-rotation coupling in changing spins. For example. For B ~ 4800 gauss, the inelastic 
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rate constant for the J Mj) = |0 1 — 1) exit channel exceeds that for the elastic one. 
This indicates that the high field limit fine structure changing collisions represent a potential 
limitation for the success of the collisional cooling processes. 



B. Collision energy dependence 

In this section, the dependence of elastic and inelastic rate constants on the collision 
energy are analyzed for the zero field limit and for three representative values of the field, 
namely 10, 200 and 4500 gauss. 

In Fig. 4. a, we report our results for collision energies 1/iK - IK, for the four mentioned 
values of the field. Elastic scattering, which is largely determined by s-waves in both incident 
and final channels, is weakly affected by the strength of the field. For spin-changing collisions, 
however, the energy dependence changes dramatically, in accordance with (|1^), which is 
shown in the figure using dashed lines. Again the trends are well-represented, although the 
formula overestimates the rates due to the zeros in the real rates described in the previous 
section. 

As a general trend, we observe that when the field is increased, the low-energy inelastic 
rate constants are substantially pushed up towards the elastic ones, but at higher energies 
the rates are less sensitive to the field. This is better illustrated in Fig. 4.b, where the 
total loss rate constants (that is the sum of the inelastic scattering in the two exothermic 
channels) for the different values of the field are plotted together with the elastic channel 
results in zero field (which, as stressed before, are essentially independent of the B value). 

In the range of collision energies from 1/iK to 1 K and for low magnetic fields the 
rates for elastic collisions remain significantly higher than the rate for the spin-flipping. In 
particular, at buffer-gas-cooling energies of ~ IK and below, the inelastic rate constants are 
in the 10"^"^ cm'^/sec range, and remain so even in the presence of a field. This result verifies 
the suitability of the ^^©2 molecule for BGC. Given the comparatively small uncertainties 
in the PES the results shown here are probably fairly realistic for the He-02 system. 
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We have continued the analysis in the range of colhsion energies from 1 to 10 K. We 
note that for energies larger than the height of the centrifugal barrier, the approximation of 
including only = 1 in the calculation is no longer accurate within a few percent, as was 
the case at lower energy. We checked that in this energy range the full calculations including 
all the possible Ai values provides results which differ (in the worst case) by about a factor of 
3 for the inelastic channels and by about a factor of 5 for the elastic one. The full calculation 
is computationally very expensive for i? > 0, as opposed to the case of zero field where a 
total- J7 basis can be adopted. Results indicate that inelastic rates for high collision energies 
are not much lower than elastic rates, and of course become higher still at energies where 
resonances exist. 

At collision energies above IK both Feshbach and shape resonances appear in the cross 
sections, as noted in p[. We find that these resonances move somewhat as a function of 
magnetic field. Nevertheless, they are sufficiently narrow that they are completely "washed 
out" by thermal averaging in a gas. We therefore present these results as a function of 
temperature rather than energy. To this end we assume a Maxwellian velocity distribution 
of the collision partners characterized by a kinetic temperature T. The thermally averaged 
rate constants are then expressed as: 



where ks is the Boltzmann constant and cr{E) stands for cross sections. To compute this 
average, the values of the cross sections for E > 10 K are extrapolated from their values at 



Averaged rate constants are shown in Fig. 5 for the same set of field values as in Fig. 4. 
The condition for magnetic trapping to be successful is usually expressed as Kei > 10 Kioss 



at least for temperatures up to 1 Kelvin and values of the field for which the asymptotic 
Zeeman splitting does not exceeds the height of the exit channel centrifugal barrier. 




(14) 



10 K. 



27| , so that we can conclude that for collisions of ^^©2 with '^He this condition is fulfilled 
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IV. APPLICATIONS 



Even in the absence of detailed information on cold collisions of a particular molecule, 
the fitting formula ([T3| ) can be used as an approximate guide to what the rates might be. 
For example, we can inquire about the prospects for evaporatively cooling ^''©2 molecules 
once they have been successfully cooled in a first stage of BGC. For this system, the d-wave 
centrifugal barrier height is ~ 13 mK, meaning that the threshold law is expected to hold 
for field values smaller than ~ 53 gauss for transitions that produce one or more molecules 
in the |0 1 — 1) final state. 

In the case of ^^©2 cold collisions we have access to the zero-field calculations of Ref. 
We have fit the energy dependence of all the inelastic rate constants to yield the dashed 
curve labeled "5 = 0" in Figure 6, which represents the total loss. The full calculation (solid 
line) has some additional features due to scattering resonances near zero energy, but this 
will not affect our conclusions here. Based on the zero-field fit, we use (0) to estimate the 
loss rates in nonzero field. The general trends are the same, namely, the rates rise sharply 
at low energy, but are roughly field-independent at larger energies. 

For evaporative cooling to be successful requires, roughly speaking, that the ratio of 
elastic to inelastic collision rates, K^i/Kioss, exceed 100 [^. For the estimated results 
shown in Figure 6, this condition holds only at energies below ~ ImK for fields as low as 
10 gauss, and not at all for near-critical fields of ~ 50 gauss. Thus evaporative cooling 
from an initially buffer-gas-cooled sample may prove trickier than previously expected. On 
the other hand, BGC is characterized by a large number of molecules cooled in the initial 
step; it is possible that a certain loss can be tolerated, and that a final sample at sub-yuK 
temperatures will still hold enough molecules to reach critical phase space density for BEG. 
Detailed kinetic simulations are required to determine if this is so. Alternatively, a recent 
proposal suggests that NH molecules could be cooled via Stark slowing to temperatures as 



low as 1 mK ||29|. In this case EG may work quite well. 

For many systems of interest to ultracold studies, there does not exist any information 
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on spin-changing collisions at low enough temperatures. In such cases it may be possible to 
make order-of-magnitude guesses anyway. For example, suppose a rate constant is known 
for higher temperature collisions. In the absence of any other information we could simply 
assert that the rate has the same value at the energy Eq corresponding to the height of the 
centrifugal barrier. The fitting formula (^) then gives the behavior of this rate at lower 
temperatures. For example, the zero-field rate constant in the first panel of Fig. 5 has a 
value of 10~^^ cm'^/sec at a temperature of 4 K. Extrapolating this value down to Eq = 0.59 
K yields a coefficient Kq = 7.2 x 10"^"^ cm^/sec for the total loss, just a factor of two higher 
than the fit to the low-energy calculation. 

V. CONCLUSIONS 

One of the main aims of this paper is to understand in a broad sense collision of param- 
agnetic molecules in the limit of ultralow temperatures in the presence of a magnetic field. 
This is a part of our effort in showing that molecules with nonzero spin in the lowest energy 
state (such as ^^©2 investigated here) can be successfully cooled and used for BEG purposes. 

Elastic and inelastic scattering in presence of a magnetic field for the specific system 
^He - ^^©2 has been characterized in detail. Our attention has been focused on the lowest 
lying trappable state of the molecule. This information is immediately relevant for BGC of 
molecular oxygen, and suggests that the presence of the field does not particularly hinder 
the BGG process. This work extends previous predictions which referred to the oxygen 
molecule in zero field |||,|^, and definitely assess the theoretical possibility for trapping this 
species. 

Moreover, we have illustrated the simple underlying physics of spin-changing rates in 
general. For fields low enough that both the incident and exit channels are in the threshold 
regime, the rate constants vary according to the Wigner-law dependence {E + AMjgfioB)^^"^ . 
This insight allows us to make estimates of rate constants beyond the ones calculated in 
detail. In particular, a small field was found to have a dramatic effect on the evaporative 
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cooling of ^ 02, which must be taken into account in future experiments aimed at quantum 
degenerate molecular gases. 
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FIGURES 

FIG. 1. The lowest-energy Zeeman levels of O2 for the even-A^ rotational manifold. These 
levels are usefully labeled by the approximate rotational (A^) and total spin (J) quantum numbers, 
along with the projection Mj of total spin onto the magnetic field. The heavy line indicates the 
lowest-lying trappable state J Mj) = |0 1 1). 

FIG. 2. Rate constants (logarithmic scale) for collisions of ^He and ^^©2 molecules in the 
J Mj) = |0 1 1) initial state as a function of the magnetic field B and collision energy 1/iK, for 
low values of the field. Solid lines are the complete quantum mechanical calculation, while dotted 
lines are results of the distorted wave Born approximation (DWBA) . For each curve the final state 
of the oxygen molecule is indicated. 

FIG. 3. Rate constants for collisions of ^He and "'^'^02 molecules in the J Mj) = |0 1 1) 
initial state as a function of the magnetic field B and collision energy 1/uK (bilogarithmic scale). 
Solid lines refer to the exact quantum calculation, while dashed lines refer to the fitting using eq. 
([1^). For each curve the final state is indicated. Vertical arrows indicate the values of the field for 
which the difference of the asymptotic energies |011) — |01 — 1) (or — |0 1 0)) exceeds the height 
of the centrifugal barrier. 

FIG. 4. (a) Collision energy dependence of elastic and inelastic rate constants in the range 1/liK 
-IK (bilogarithmic scale). In each panel the corresponding value of the magnetic field is indicated. 
The solid lines refer to the full quantum calculations while dashed lines refer to the inelastic rate 



constants calculated using relation (13), with values of Kq as given in the text, (b) Shows the 
total loss rate constant for all four field values simultaneously, to facilitate their comparison. In 
(b), only the full quantum calculations are shown. 

FIG. 5. Thermally averaged elastic (solid line) and total loss (dashed line) rates for ^He — ^^02 
collisions as a function of temperature. In each panel, the corresponding value of the magnetic 
field is indicated. 
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FIG. 6. Comparison between elastic and total loss rate constants for collisions of O2 molecules 
with each other, for several representative values of the field. The zero-field results (solid line) are 
reproduced from Ref. [9]. The field-dependent estimates (dashed lines) are based on a fit of eq. 
( [T^ to the zero-field rates. Fitted results are shown only for field values and collision energy ranges 
where eq. ( p^ ) is expected to hold. 
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